N89 - 19260 


AIRFOIL STALL PENETRATION AT CONSTANT PITCH RATE 
AND HIGH REYNOLDS NUMBER 


Peter F. Lorber and Franklin 0. Carta 
United Technologies Research Center 
East Hartford, Conn. 


PRECEDING PAGE BLANK NOT FILMED 


519 



MODEL WING IN WIND TUNNEL 


The model wing consists of a set of fiberglass panels mounted on a steel spar that 
spans the 8 ft. test section of the UTRC Large Subsonic Wind Tunnel. The pitch angle 
of the model is set by hydraulic actuators attached to each end of the spar and 
controlled using a dual channel closed-loop system. This allows 1024 samples of each 
data channel to be acquired at identical positions during each cycle of the wing 
motion. The first use of this system was to measure surface pressures and flow 
conditions for a series of constant pitch rate ramps and sinusoidal oscillations at 
Mach numbers between 0.2 and 0.4, Reynolds numbers between 2 and 4x10^, and pitch 
rates between A = cd/2U = 0.001 and 0.02. 
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AIRFOIL SECTION 


A Sikorsky SSC-A09 airfoil with a 17.3 in. chord was used. This is a 9 % thick 
supercritical section designed for low drag at high subsonic Mach numbers. 
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INSTRUMENTATION 


The wing-mounted instrumentation consists of 72 miniature pressure transducers and 8 
surface hot film gages. The pressure transducers are mounted inside the fiberglass 
skin and are connected to the surface by short pipettes. This technique provides a 
point measurement on a smooth surface contour. The primary chordwise pressure 
transducers are located in arrays of 18 on each surface, 0.5 chordlengths from the 
tunnel centerline, and 2.3 chordlengths from the side wall. The transducers are 
arranged in a segmented Gaussian array from x/c = 0.005 to x/c = 0.99. The other 36 
transducers are located in additional arrays at 1.8 and 1.4 chordlengths from the 
wall, intended for use in future swept wing and finite tip experiments. The hot film 
gages are located in an array parallel to the primary pressure array, from x/c = 

0.026 to x/c = 0.88. Detailed calibrations over both pressure and temperature ranges 
make it possible to measure both mean and unsteady pressures to within 1 % of the 
calibration range. An interactive data system is used to acquire single samples of 
time histories, or to obtain an ensemble average of the time histories, based on 20 
cycles of the motion. These data are digitized, converted to pressure coefficient, 
and stored on magnetic tape. 
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TEST CONDITIONS 


Unsteady data were acquired for 36 constant pitch rate ramps and 9 sinusoidal 
oscillations, as shown in the table. The maximum pitch rate of 360 deg/sec occurred 
at A * 0.02 and M * 0.2, and the minimum of 18 deg/sec at A * 0.001. The maximum 
rate is lower than the maximum reached in many smaller scale experiments (Refs. 

1-5), but is larger than that for both the "typical" (A = 0.001) and "minimum time" 
(A - 0.0044) maneuvers described in Ref. 6. The wing angle of attack was limited by 
the support system to a maximum of 30 deg, which was less than the maximum obtained 
in the smaller scale tests, but sufficient to include all of the primary 
stall-related events at the pitch rates used. 
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STEADY PRESSURES 


The steady pressure distributions have a strong leading edge suction peak prior to 
stall. Maximum suction occurs at a = 13 deg and reaches values of Cp = -8.2 at 
M * 0.2, and C = -5.2 at M =0.4. The sonic pressure coefficient at M = 0.4 is 
-3.7, so a small supersonic zone exists at the leading edge, terminated in a shock 
near x/c - 0.03. The compressibility effects associated with this zone create the 
lift curve differences seen earlier. 
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RAMP TIME HISTORIES 


Airfoil motion for a ramp consists of an initial delay 
minimum angle, a constant rate increase to the maximum 
maximum angle, and a return to the initial condition, 
a small portion, T, of this cycle. Time histories of 
acquisition period are shown in the figure for several 
The ramp begins at nondimensional time, r = t/T, of 0. 
pitch increase is quite linear with time and has sharp 
higher pitch rates the damping of the hydraulic system 
maintains a nearly linear pitch rate near stall. 


of several seconds at the 
angle, a second delay at the 
Data are acquired only during 
the pitch angle during the data 
0 to 30 deg ramps at M = 0.2. 
125 and ends at t = 0.625. The 
corners for A < 0.005. At 
rounds the corners, but still 
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UNSTEADY PRESSURES, UPPER SURFACE 


Ensemble averaged upper surface pressure time histories for 11 of the 18 measurement 
stations are shown for a 0 to 30 deg constant pitch rate ramp at A = 0.005, M = 0.2. 
The ordinate scale on the lower left refers absolutely to the x/c - 0.005 curve, and 
all other curves are vertically offset by AC - 1.0. Each of these curves is 
referenced to its own origin (upper left scale). The pressure responds smoothly to 
the imposed pitch angle until maximum suction is reached at r = 0.4 (a = 15.5 deg) at 
letter A. Peak suction at x/c - 0.005 is Cp - -12.5, corresponding to a local Mach 

number of 0.84. The rapid local increase in suction associated with passage of the 
vortex (letter B) is only observed for x/c > 0.1,, implying that the vortex forms 
ahead of this position and travels downstream thereafter. The vortex pressure 
propagates along the chord at 0.16U, and reaches the trailing edge at r * 0.45. For 
t> 0.50 the upper surface separation is massive and the pressures are virtually 
invariant in both time and position (letter C). All events are completed well before 
the ramp ends at r = 0.625. 
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UNSTEADY PRESSURES, LOWER SURFACE 


The lower surface does not separate, and responds primarily to the migration of the 
stagnation point near the leading edge (from motion initiation at t = 0.125 to r* 
0.42), and to the upper surface massive stall near the trailing edge (letter A), 
where the pressure is required by continuity to match that on the upper surface. In 
addition, a periodic oscillation in these ensemble averaged pressures is detected 
near the lower surface trailing edge (letter B). The oscillation frequency of 62Hz 
is very close to the 65Hz frequency calculated for the vonKarman vortex street shed 
by a circular cylinder having a diameter, D, equal to the vertical projection of the 
airfoil chord, c sin a. From Ref. 7, a cylinder will generate a vortex street at a 
Strouhal number, Df/U, of 0.21 for Reynolds Numbers less than 10^ and at fD/U - 0.27 
for Re > 3 x 10^. No regular vortex street is formed when the Reynolds number is 
between these limits. In the present case the Reynolds number of the cylinder is 
equivalent in projected area to the airfoil at M - 0.2 and a - 30 deg is 10^, near 
the boundary where periodic oscillations should cease. 
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INDIVIDUAL RECORDS, LOWER SURFACE 


A larger post-stall oscillation amplitude is found for individual records of the 
motion than for the ensemble average. The top time history in this figure is the 
ensemble averaged trailing edge plot from the previous figure, and the remaining 
curves are a set of (nonconsecutive) individual records at the same location and for 
the same set of flow conditions. Each curve is offset vertically by Cp — 0.5. The 
oscillations have the same frequency and a similar maximum amplitude during each 
record, but they are not well correlated in phase. The resulting cancellation 
reduces the ensemble averaged amplitude. A separate oscillation at a lower 
frequency is often present immediately following the passage of the stall vortex at 
t - 0.45 (letter A), and may result from secondary stall vortices. 
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HOT FILM TIME HISTORIES 


Hot film gages provide additional information on the surface flow conditions, as 
shown by the six time histories in this figure. The data are self-scaled, so each 
varies over one unit. As indicated on the leading edge trace, the boundary layer is 
initially laminar at the first two stations, goes through transition, becomes 
turbulent, and eventually becomes separated. From x/c = 0.192 and aft the flow is 
always either turbulent, or separated, and is never laminar. Ahead of the stall 
vortex release point (x/c * 0.1) the heat transfer drops when the boundary layer 
separates at r = 0.4. Downstream of the release point the heat transfer drops 
slightly at separation but then rises rapidly as the high velocities induced by the 
stall vortex pass each gage. There is excellent agreement between the times of 
separation as measured by the hot film gages and by the adjacent pressure 
transducers. 
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UNSTEADY PRESSURE, UPPER SURFACE, MAX PITCH RATE 


At the maximum pitch rate of A = 0.02 the upper surface pressure time histories are 
qualitatively similar to the lower pitch rate results prior to stall, but are 
significantly different after stall. The unsteady increments to the airloads are 
strongly influenced by the strength of the stall vortex, which may be estimated by 
measuring the local change in Cp at some point along the chord downstream of the 
vortex release point, say at x/c - 0.302. This increment increases from AC - 1.0 
at A - 0.001 (not shown here) to AC - 1.5 at A - 0.005 and to AC - 2.5 at P A - 
0.020 (below). p p 
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VORTEX PROPAGATION TRACES 


'The vortex propagation speed along the airfoil may be estimated using the times of 
minimum pressure at each transducer. These points are plotted below for 5 pitch 
rates. The region of constant vortex speed generally begins near x/c = 0.10 and ends 
between x/c = 0.80 and 0.90, and the speed in this region increases approximately 
linearly with pitch rates, from 0.13U at A = 0.001 to 0.33U at A = 0.020. These 
results are consistent with those of previous investigations (Ref. 8 for sinusoidal 
motions) . 
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UNSTEADY FORCE AND MOMENT COEFFICIENTS 


Integrated lift, pressure drag, and pitching moment (parts a, b, c) exhibit 
characteristic behavior with increasing pitch rate: (1) Before the stall vortex is 
formed, increasing pitch rate decreases the lift slope, decreases pitching moment, 
and increases drag. The lift slope effect agrees with the results of Ref. 2 at lower 
Reynolds number. (2) There is a rapid buildup of lift as the leading edge vortex 
forms. The unsteady increment added to the quasi-steady lift increases from AC^ * 

0.4 at A = 0.001 to 1.1 at A = 0.020. (2) The pressure drag increases smoothly 

before stall, rises and falls rapidly as the stall vortex travels over the chord, 
then increases slowly as the pitching of the airfoil rotates the aerodynamic force 
vector. (4) The pitching moment in attached flow becomes more negative at higher 
pitch rate, following the prediction of thin airfoil theory (cf. Eq. 4-171 in Ref. 

9). (5) The peak negative C M after stall increases with pitch rate, from C M = 

0.22 at A = 0.001 to -0.52 at A > 0.020. 
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SEQUENCE OF STALL EVENTS 


The sequence of stall events is systematic and each appears to be linear with pitch 
rate. First the stall vortex forms and and start to rise. Moment stall 
occurs when the vortex is released. Maximum occurs as the vortex travels 
downstream along the chord, and when the vortex reaches the trailing edge, the 
minimum and maximum C D are obtained. At these pitch rates the angle when moment 
stall occurs increases approximately linearly with A, in agreement with results for 
a sinusoid at similar peak pitch rates (Ref. 8). This does not agree with the 
square root correlation postulated by Gormont (Ref. 10). 
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UNSTEADY PRESSURES, UPPER SURFACE, MAX MACH NO. 


A supersonic zone at the leading edge has a profound effect on chordwise pressures at 

freestream Mach number M = 0.4. The time histories for this 0 to 20 deg ramp at A - 

0.005 are not shifted vertically and the ordinate is an absolute measure for all 
curves. At x/c = 0.005 the flow becomes supersonic at r = 0.34 (a = 7.6 deg) and 
reaches a peak local Mach number of 1.3 at r = 0.45 (a* 10.6 deg). The sharp rise 
and fall of the pressure at x/c = 0.026 (letter A) is caused by the movement of the 
shock downstream past this station at r = 0.40 and back upstream at r = 0.46. A 
distinct stall vortex is released at r = 0.46 and x/c = 0.060, just downstream of the 
shock (letter B). Several of the vortex characteristics differ from those at M • 

0.2: (1) the release point is at x/c = 0.060 instead of 0.10; (2) the vortex speed 

is approximately 10% less than the speed at M = 0.2; (3) the strength of the vortex 
is reduced by approximately 50%; and (4) the pressure signature of the vortex is not 
observed downstream of x/c = 0.57. In addition, no clearly defined oscillations are 
present after stall at M - 0.4. The Reynolds number of the equivalent bluff body is 

1.4 x 10^ for this case, above the maximum of 1 x 10^ for a stable vortex street 

(Ref. 7). 
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Lift Force, 







UNSTEADY PRESSURES, UPPER SURFACE, SINUSOID 
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Sinusoidal motion differs from ramp motion in two important aspects: in a sinusoid 

the airfoil does not start from a steady-state condition, and the pitch rate is 
constantly changing. Nonetheless, there are many qualitative similarities between 
the responses to the two motions. For the upper surface time histories shown here 
(a = 20 - 10 coso»t at M = 0.2 and k « 0.05) the motion begins at r = 0.0, reaches 
maximum pitch angle at t « 0.50, and returns to the minimum value at r * 1.0. Each 
time history here is offset vertically by AC = 1.0. The characteristics are 
similar to the constant pitch rate results sRown earlier for A = 0.005, which is 
comparable to the instantaneous sinusoidal pitch rate at stall of A = 0.008. 
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UNSTEADY PRESSURES, LOWER SURFACE, SINUSOID 


The lower surface pressures for this same sinusoid are remarkable in that the maximum 
pressure change occurs at the trailing edge (to match the stalled suction surface 
pressure level), exhibits an initially damped behavior, and then stablilizes into a 
coherent vortex street response. The latter is noteworthy because in this 
ensemble-averaged data set the vortex street frequency must be synchronized with the 
fundamental frequency of the sinusoid. No single records were taken to permit close 
examination of the vortex street oscillation. 
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TRANSITION POINT DETECTION 


Accurate determination of the point of boundary- layer transition is critical to many 
numerical computations. Surface heat transfer gages are usually required to 
determine the unsteady motion of the transition point. These gages are generally 
cumbersome to install and each requires a dedicated active anemometer circuit to 
operate. The techniques developed during the course of this study show that it may 
be possible to locate the transition point using the local increase in rms pressure 
measured by the unsteady pressure transducer as the local flow changes its character. 
An example of the correlation between the hot film and pressure transducer results 
is shown below. All results are at x/c « 0.026 for a 0 to 30 deg ramp at A - 0.02 
and M = 0.20. The vertical pressure scale has been severely foreshortened, and only 
a small portion of the time scale is displayed. The increase in pressure rms 
corresponds to the increase in ensemble averaged heat transfer (upper curve) and 
to the spike in the rms heat transfer (second curve). No change at transition can 
be discerned in the ensemble-averaged pressure (bottom curve). If these promising 
early results are confirmed by additional correlations, this method may make it 
easier to locate transition in complex three-dimensional and unsteady flows. 
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NOMENCLATURE 


A pitch rate, ac/2U, rad/sec 

c airfoil chord, m 

Cp section pressure drag coefficient, D/Qc 

C. section lift coefficient, L/Qc 

L 2 
section pitching moment coefficient about x/c * 0.25, m/Qc 

Cp pressure coefficient, (P-Pstatic)/Q 

k reduced frequency of sinusoidal motion, o»c/2U 

2 

Q freestream dynamic pressure, 0.5 p U , Pa 

t time from start of data acquisition, sec 

T data acquisition period, sec 

Re Reynolds number, c U/i> 

U freestream velocity, m/sec 

x distance from airfoil leading edge, m 

a geometric pitch angle, deg 

2 

v kinematic viscosity, m /sec 

3 

p air density, kg/m 

r nondimensional time, t/T 

a> frequency of sinusoidal oscillation, 2 n f 


540 



CONCLUSIONS 


Conclusions - Increased Pitch Rate 

• Stall events are delayed 
•Stall vortex is strengthened 
•Vortex propagation speed increases 
•Unsteady airloads increase 


Conclusions - M = 0.3, 0.4 

• Supersonic zone near leading edge 

• Stall vortex is weaker 

• Unsteady airloads are reduced 

Additional Conclusion 

• Post stall vortex shedding when Re bi u ff < 10® 

• Shedding synchronized when A > 0.01 

• Sinusoids and ramps qualitatively similar 


Note: This work was supported by the U.S. Air Force Office 

of Scientific Research under Contract F49620-84-0082 
and will be reported in Ref. 11. 
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